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Abstract 

The macroscopic control of ubiquitous heat flow remains poorly explored due to the lack of 
a fundamental theoretical method. Here, by establishing temperature-dependent transformation 
thermotics for treating materials whose conductivity depends on temperature, we show analytical 
and simulation evidence for switchable thermal cloaking and a macroscopic thermal diode based 
on the cloaking. The latter allows heat flow in one direction but prohibits the flow in the opposite 
direction, which is also confirmed by our experiments. Our results suggest that the temperature- 
dependent transformation thermotics could be a fundamental theoretical method for achieving 
macroscopic heat rectification, and provide guidance both for macroscopic control of heat flow and 
for the design of the counterparts of switchable thermal cloaks or macroscopic thermal diodes in 
other fields like seismology, acoustics, electromagnetics, or matter waves. 
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Heat flow is a ubiquitous phenomenon in nature, and hence how to control the flow 
of heat at will is of particular importance for human life. Fortunately the past years have 
witnessed the possibility of manipulating phononic [l-7| or electronic js-lll heat conduction 
at the nanoscale, which provides a promising method to make use of heat flux. So far 


signihcant progresses 


and caloritronics js, 9|. On the contrary, steering heat conduction at the macroscale is 


rave been made, such as computation j^, information storage jd] 


still far from being satisfactory 


12| . which prohibits macroscopic thermal rectihcation for 


diverse thermal management problems, e.g., efficient refrigerators, solar cells, and energy¬ 


saving buildings 


13l-ll5|. A reason behind the situation is due to the lack of a fundamental 


theoretical method. 

In 2008, Fan et al. adopted a coordinate transformation approach to propose a class 
of thermal metamaterial where heat is caused to flow around an “invisible” region at steady 
state, thus called thermal cloaking. The cloaking originates from the fact that the thermal 
conduction equation remains form-invariant under coordinate transformation. So far, the 


theoretical proposal of steady-state thermal cloaking m and its extensions fsay, bifunctional 
cloaking of heat and electricity 17|] or nonsteady-state thermal cloaking ISj) have been ex- 

19H23| . The theoretical treatment based on coordinate 


perimentally verihed and c 
transformation 


16-18 


24 


eveloped 


28| . which is called transformation thermotics (or transformation 


thermodynamics) , has a potential to become a fundamental theoretical method for macro- 
scopically manipulating heat flow at will. 

However, in order to realize desired macroscopic thermal rectihcation, say, switch- 
able thermal cloaking and macroscopic thermal diodes, the existing transformation ther¬ 
motics 


16 


18[ | is not enough since it only holds for materials whose conductivity is inde¬ 


pendent of temperature (thus called linear materials). For instance, the desired thermal 
diode should conduct heat in one direction, but insulate the heat in the opposite direction. 
This is a kind of asymmetric behavior of heat current. For this purpose, nonlinear materials 
(whose conductivity relies on temperature) must be adopted. Actually, it is long known 
that for many materials, their thermal conductivities {k) vary with temperature (T): (1) 
K increases as T increases. Say, for noncrystalline solids, a series of experiments on glass 
29l | showed that at low temperature the thermal conductivity is proportional to ® ~ 


(2) K increases as T decreases. For example, measurements on single crystals of silicon 
and germanium from 3 K to their melting point 30j showed that their thermal conductivity 
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decreases faster than 1/T. 

Temperature-dependent transformation thermotics 

Now we are in a position to establish a theory of transformation thermotics for treating 
nonlinear materials, thns called temperature-de pen dent transformation thermotics. The de¬ 
tails of the theory are given in Part I of Ref. 3^, which yield the key formula, Eq. (S8). 
This Eq. (S8) denotes that instead of constructing materials for a background whose ther¬ 
mal conductivity is T-dependent, we may apply a T-involved transformation to the original 
thermal conduction equation. This process allows us to design switchable thermal cloaks 


and then macroscopic thermal diodes. The 
investigated cloaks without switches 


16 


231, 


brmer serve as an extension of the extensively 


3l| : the latter are actually a useful application 


of the former in this work. We proceed as follows. 

Switchable thermal cloaks 

Design .— A traditional thermal cloak can protect a central region from an external 
heat flux and permits the region to remain a constant temperature without disturbing the 
temperature distribution outside the cloak (Fig. 1). To achieve this cloaking effect, a simple 
radial stretch transformation of polar coordinates may be performed. As schematically 
shown in Fig. 1, the circular region with radius i ?2 is compressed to the annulus region with 
radius in-between Ri and i? 2 , and the geometrical transformation can be written as 

i?2 ~ Ri 




Ri 


+ Ri: 


( 1 ) 


where r G [0, Rfl and r' G [i?i, Rfl. Here r' is the radial coordinate in physical space. 

In order to realize different responses to heat flow on the two sides of a thermal diode, we 
need two types of thermal cloaks: one functions at high temperature (hereafter indicated as 
type-A cloaks), the other works at low temperature (type-B cloaks). Unlike some previous 


proposed switchable electromagnetic cloaks j32l-l34| . the switching effect should be triggered 


automatically by temperature changes. For this purpose, an idea is to modify Eq. o as 

Ri-Ri{T) 


y .' — ^"^2 


Ri 


+ Ri(T), 


( 2 ) 


where Ri{T) = i?i[l — (1-|-for type-A cloaks and .Ri(T) = i?i/(l-1-for 
type-B cloaks. Here is a critical temperature around which the cloak is switched on or 
off, and /9 is a scaling coefficient which is set to be 2.5 in this work. 
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So far, for obtaining thermal cloaks with switching phenomena, we need to combine 
Eq. ([2]) and Eq. (S8). As a result, for the area with radius r' e [i?i, R 2 ] in Fig. 1, we achieve 
the thermal conductivities in polar coordinates, diag[kr(T),ke(T)], as 


Kr{T) = Kq- 


R,{T) 


, K0{T) = Ko- 


, 4-. (3) 

r r'-Ri{T) 

Here Kq represents the T-independent thermal conductivity of the background. 

Finite-element simulations .— Then we perform hnite element simulations based on the 
commercial software COMSOL Multiphysics. A thermal cloak with i?i = 1 cm and R 2 = 
2 cm is set in a box with size 8 cm x 7 cm as shown in Fig. 2. Heat diffuses from the 
left boundary with high temperature Th to the right boundary with low temperature T^. 
Meanwhile, the upper and lower boundaries of the simulation box are thermally isolated. 

The simulation results of a type-A cloak are shown in Fig. 2(a,b). In Fig. 2(a), at 
high temperature (T = 340 K ~ 360 K), we observe that the cloak is functioning and 
thermally hiding the object located at the central region with radius Ri. However, when 
the environment is changed to low temperature (T = 300 K ~ 320 K), the cloak is turned 
off. That is, the temperature distribution outside the object is distorted, just as the cloak 
(located between Ri and R 2 ) is absent. Owing to the antisymmetry between type-A and 
type-B cloaks, the type-B cloaks exhibit the behavior similar to Fig. 2(a,b), but switching 
on (or off) at low (or high) temperature. 

Theoretical realization based on the effective medium theory .— The materials designed 
according to Eq. ([3]) is anisotropic and inhomogeneous, which is difficult to be realized in 
experiments. In fact, for constructing such materials, we can simply utilize alternating layers 
of two homogeneous isotropic sub-layers with thicknesses di and ^2 and conductivities Ka and 
Kb. For simplihcation, in this work, we set d] = (R. According to the theoretical analysis 


and effective medium theory (EMT) 


17 


19 


31 


35|, the conductivities of two sub-layers 


should satisfy KaKb ~ k^ for a traditional cloak. In order to endue conventional cloaks with 
the switching effect, some mathematical operations must be carried out on Ka and Kb in the 
way analogous to what we did on Ri{T). Therefore, we obtain ki{T) and K 2 {T) as the new 
temperature-dependent thermal conductivities of the sub-layers. 


^ QpiT-Tc) ’ I qRT-Tc) ' 

The two expressions yield: ki{T) — )■ Ka and K2(T) — )■ Kb when T ^ Tc; ki(T) kq and 
K 2 iT) kq as T -C Tc- (The detailed dependence of Ki{T) and K 2 (T) on temperature 
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is displayed in Fig. S3 of Part III of Ref. [3^.) That is, the cloaking effect is switched on 
(or off) for high (or low) temperature environment T ^ Tc (or T <C Tc). Eq. (j3]) offers a 
convenient tool to help experimentally realize our theoretical design of Eq. ([3]). Next, we plot 
Fig. 2(c,d), which shows the simulation results of 10 alternating layers for constructing type- 
A cloaks. Evidently, Fig. 2(c,d) displays the switching phenomenon similar to Fig. 2(a,b). 
The procedure holds the same for achieving type-B cloaks. 

Macro thermal diode 

Design .— The above thermal cloaking may help to design a kind of macroscopic thermal 
diodes. As shown in Fig. 3(a), the diode device contains Regions I, II and III; Region I (II) 
is a segment of the type-A (type-B) cloak, and Region III is a thermal conductor. Compared 
with a full system, the cloaking effect still exists in our diode but is not perfect. There will 
be a small amount of heat flux conducting through the central region for the insulating case. 
However, the truncation of a whole cloak is necessary to separate the type-A part and the 
type-B part. The antisymmetry of type-A and type-B cloaks is expected to cause different 
behaviors of heat conducting from the two opposite directions. The transformation plays an 
important role in introducing anisotropic effect to the structure, which guides the heat flux 
to the boundaries to enhance the thermal insulating effect. 

Finite-element simulations .— Then we perform finite element simulations. Fig. 3(c,d) 
shows the simulation results of the device, which helps to insulate heat from right to left 
but conduct the heat from left to right. That is, the behavior of diode can be achieved 
indeed due to the antisymmetry of type-A and type-B cloaks (namely. Region I and Region 
II). Moreover, for different temperature biases (obtained by subtracting the temperature at 
the left boundary from that at the right boundary), AT, we also calculate the total heat 
current J by integrating the x component of heat flux across the line x = 0; see Fig. 3(b). 
The device displays a significant rectifying effect, which has a maximum rectification ratio 
of 30 for the current parameter set. 

Experimental realization based on the effective medium theory .— In order to realize such a 
macroscopic thermal diode, we can also adopt the EMT. As discussed above, the switchable 
thermal cloak can be constructed with alternating layers. The thermal conductivities of the 
sub-layers are required to be sensitive to the temperature change around the critical poin 


This kind of behaviors can be found in the phase transitions of some materials 


37 


39| 


However, inspired by the spirit of metamaterials (yielding novel functions by assembling 
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conventional materials into specific structures), we manage to realize the macroscopic ther¬ 
mal diode with materials of constant conductivities. Our method is that instead of directly 
resorting to the transitions of materials’ physical properties, we use the structural transition 
to trigger the switching effect. According to the demands of our design, the geometrical 
conhguration of the device should change rapidly as temperature varies. Fortunately, we 
found that the shape-memory alloy (SMA) ^ may be able to help. As shown in the 
schematic diagrams of our desi gn [see Fig. 4(a) and (b); more details of our experiment 


can be found in Part II of Ref. 


36l|]. the sub-layers of the cloak segments are coppers and 


expanded-polystyrene (EPS). Around the critical temperature, the deformations of the SMA 
slices drive the copper slices to connect or disconnect the copper layers. The connection and 
disconnection can be equivalently regarded as transitions of the local thermal conductivities. 
Thus a temperature-dependent thermal metamaterial is realized with materials of constant 
thermal conductivities (a whole switchable thermal cloak can also be built with the same 
method). The experimental results are shown in Fig. 4(c) and (d). Fig. 4(c) displays the 
temperature distribution within the central region, which is almost constant. That is, in this 
case, heat almost cannot flow through this central region. Thus, Fig. 4(c) corresponds to 
the insulating case of the diode. In contrast. Fig. 4(d) represents the conducting case of the 
diode. This is because Fig. 4(d) shows a significant temperature gradient within the central 
region. Also, the temperature distribution appears to be horizontal. As a result, we can 
conclude that an evident heat flow comes to appear within the central region of Fig. 4(d). 

Conclusions 

We have established a theory of temperature-dependent transformation thermotics for 
dealing with thermal materials whose conductivity is temperature-dependent. The theory 
serves as a fundamental theoretical method for designing switchable thermal cloaking. We 
have also shown that the switchable thermal cloaks can be employed for achieving a macro¬ 
scopic thermal diode, which has also been experimentally realized by assembling homoge¬ 
neous and isotropic materials according to the design based on the EMT (effective medium 
theory). The diode has plenty o^otential applications related to heat preservation, heat 


dissipation, or even heat illusion 


42| in many areas like efficient refrigerators, solar cells. 


energy-saving buildings, and military camouflage. Thus, by using temperature-dependent 
transformation thermotics to tailor nonlinear effects appropriately, it becomes possible to 
achieve desired thermal metamaterials with diverse capacities for macroscopic thermal rec- 
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tification. On the same footing, our consideration (for cloaks and diodes) adopted in this 
work can be extended to obtain the counterparts of both switchable thermal cloaks and 
macroscopic thermal diodes in other helds like seismology, acoustics, electromagnetics, or 
matter waves. 
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Fig. 1 (color online). Schematic graph depicting a thermal cloak between radius Ri and 
/? 2 - Red lines with arrows denote the flow of heat: the cloak does not disturb the heat flow 
at the region with radius larger than i? 2 ; the heat flux cannot enter the central region with 
radius smaller than Ri. 

Fig. 2 (color online). Switchable thermal cloaks obtained by two-dimensional finite- 
element simulations: (a,c) switch on for the temperature above 340 K and (b,d) switch off 
for the temperature below 320 K. The color surface denotes the distribution of temperature, 
where isothermal lines are indicated; heat diffuses from left to right; the upper and lower 
boundaries are thermal insulation, (a) and (b) show the results of thermal conductivities 
calculated according to Eq. ([2]); (c) and (d) show the results of 10 alternating layers of 
two sub-layers with Hi{T) and H 2 iT) given by Eq. (jl]) (EMT). In (a-d), an object with 
thermal conductivity 0.01 W/mK is set in the central region with radius Ri. Parameters: 
kq = 1 W/mK, Kq = 0.1 W/mK, ki, = 10 W/mK, Ri = 1 cm, R 2 = 2 cm, and Tc = 330 K. 

Fig. 3 (color online), (a) Sketch of a thermal diode, which is the rectangular area enclosed 
by the solid black lines. The blurred area outside is a reference and actually does not 
exist in the design. I, II, and III represent three regions, respectively. Here the arrows 
indicate the direction of heat flow; the length of arrows represents the amount of heat 
flux: the heat flux transferred from right to left (upper panel: the insulating case) is much 
smaller than that from left to right (lower panel: the conducting case), (b) Heat current 
J versus temperature bias AT. (c,d) Thermal diode obtained by two-dimensional finite- 
element simulations: (c) the insulating case and (d) the conducting case. The color surface 
denotes the distribution of temperature; white arrows represent the direction of heat flow; 
the length of white arrows indicates the amount of heat flux; the upper and lower boundaries 
are thermal insulation. Thermal conductivities are calculated according to Eq. (|3]); an object 
with thermal conductivity 10 W/mK is set in the central region with radius Ri. Parameters: 
Ko = 1 W/mK, Ri = 3.6 cm, R 2 = 4 cm, and Tc = 330 K. 

Fig. 4 (color online). (a,b) Scheme of experimental demonstration of the macroscopic 
thermal diode: (a) insulating case and (b) conducting case. Both the copper-made concentric 
layered structure and the central copper plate (both displayed in orange) are placed on an 
expanded-polystyrene (EPS) plate which is not shown for clarity. The left and right sides 
of the diode are stuck in water to have constant temperature boundary conditions, (a) 
When cold water is filled in the left container (light blue) and hot water the right container 
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(pale red), the bimetallic strips of SMA and copper (white) warp up and the device blocks 
heat from right to left, (b) When the two containers swap their locations, the bimetallic 
strips (white) flatten and the device conducts heat from left to right. (c,d) Experimentally 
measured temperature distribution of the device: (c) insulating case and (d) conducting 
case. 
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FIG. 1: /Li, Shell, Wu, Huang, Chen, Ni, and Huang 
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FIG. 3: /Li, Shen, Wu, Huang, Chen, Ni, and Huang 
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FIG. 4: /Li, Shen, Wu, Huang, Chen, Ni, and Huang 
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